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1.  INTRODUCTION 

The  purpose  of  this  paper  is  to  propose  and  investigate  a  flow  con¬ 
trol  algorithm  for  adjusting  session  rates  in  a  data  or  voice  network. 

This  algorithm  is  motivated  by  a  Voice  Coder  scheme  introduced  in  [1] 
and  called  "embedded  coding".  In  this  scheme  a  segment  of  talkspurt 
is  coded  into  packets  of  different  priority  levels.  The  higher  priority 
packets  contain  the  "core"  of  the  speech  while  the  lower  priority  packets 
contain  the  information  that  "fine  tunes"  it.  Traditional  voice  flow 
control  mechanisms  either  block  the  initiation  of  a  call  or  discard  small 
segments  of  it  while  it  is  already  in  progress.  By  contrast  the  embedded 
coding  scheme  dynamically  trades  off  between  voice  quality  and  congestion 
by  discarding  the  lower  "priority"  packets  either  at  the  point  of  congestion 
or  the  point  of  entry.  The  level  of  congestion  at  which  the  gaps  between 
the  segments,  delivered  by  the  traditional  schemes,  render  the  speech  un¬ 
intelligible  is  much  lower  than  the  one  at  which  the  embedded  coding 
scheme  delivers  unintelligible  information.  This  flexibility  in  exercis¬ 
ing  flow  control  makes  the  embedded  coding  scheme  attractive. 

Alleviation  and  prevention  of  congestion  by  discarding  lower  priority 
packets  at  the  point  of  entry  seems  to  be  superior  to  discarding  them  at 
1  the  point  of  congestion.  The  latter  amounts  to  a  waste  of  network  resources. 

|  But,  it  would  not  be  advisable  to  forgo  the  capability  of  discarding  lower 

i 

j  priority  packets  at  the  point  of  congestion,  because  of  the  time  delay 

|  involved  in  making  the  entry  points  aware  of  downstream  congestion.  Based 

j  on  this  we  believe  that  both  capabilities  should  be  used.  The  rates  at 
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the  entry  points  will  be  determined  upon  longer  time  averages  of  congestion 
levels  while  the  capability  of  discarding  packets  at  the  point  of  congestion 
will  serve  to  alleviate  intolerable  momentary  bottlenecks.  The  rates  at 
the  entry  point  will  be  adjusted  so  that  the  capability  of  discarding 
packets  at  the  point  of  congestion  will  not  be  exercised  too  often. 

In  this  paper  we  discuss  a  method  of  determining  the  input  rates  at 
the  entry  point.  To  this  end  we  will  ignore  the  capability  of  discarding 
packets  within  the  network  in  order  to  simplify  the  analysis.  As  in 
quasistatic  routing  we  employ  an  "on-line"  iterative  algorithm  that  will 
solve  a  static  problem.  The  hope  is  that  the  algorithm  converges  fast 
enough  relative  to  the  session  initiation  and  termination  process,  and 
as  a  result  will  be  able  to  "track"  its  variation  keeping  the  rates  in 
the  ballpark  of  the  optimal  rates  at  all  times. 

The  criterion  used  to  determine  input  rates  is  based  on  the  notion 
of  "fair  allocation"  introduced  in  Section  2.  Roughly  speaking  the 
objective  is  to  maximize  the  smallest  session  rate,  and  once  this  is 
achieved  to  maximize  the  second  smallest  rate,  etc.  In  Section  3  we 
introduce  the  algorithm  and  describe  its  convergence  properties. 

The  idea  of  the  algorithm  is  to  adjust  the  input  rates  of  sessions 
on  the  basis  of  the  current  level  of  congestion  along  the  session  path. 

The  necessary  information  is  collected  by  a  control  packet  sent  period¬ 
ically  by  each  session  origin  along  the  path  similarly  as  in  flow  control 
methods  investigated  by  simulation  in  [1].  This  method  of  adjusting 
input  rates  seems  also  suitable  for  other  situations  where  fast  reaction 
to  momentary  congestion  is  needed.  For  example  when  the  number  of  users 


ft 


in  the  network  is  small  but  some  of  these  users  can  overload  the  network 
if  left  uncontrolled  to  transmit  at  maximum  rate,  then  a  dynamic  method 
of  routing  and  flow  control  is  needed.  The  ideas  of  this  paper  can 
provide  an  alternative  to  existing  techniques  [7],  [8]  in  such  situations. 
Further  work  is  required  along  this  direction. 


I 
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2.  PROBLEM  FORMULATION 

Consider  a  network  with  nodes  1,  2,  . ..,  N  and  a  set  of  directed 
links  L •  Each  link  a cL  has  a  capacity  C&  associated  with  it--a  positive 
number.  Let  S  denote  a  set  of  sessions  taking  place  between  nodes.  Each 
session  seS  has  an  origin  node  associated  with  it  and  traverses  a  subset 
of  links  denoted  by  Lg.  Note  that  we  do  not  restrict  the  session  to  have 
a  single  destination,  so  the  set  of  links  Ls  may  be  for  example  a  tree  rooted 
at  the  origin  node  of  s  and  used  for  broadcasting  messages  throughout 
the  network.  We  denote  by  S  the  set  of  sessions  traversing  a 
link  aeL.  If  rg  is  the  input  rate  of  session  s  (in  data  units/sec),  then 
the  flow  F  of  a  link  aeL  is  given  by 

cl 


F 

a 


l 

Se^a 


V 


(1) 


The  problem  broadly  stated  is  to  choose  a  vector  of  session  input  rates 
r  =  (...,  rs>  ...)  which  results  in  a  set  of  "satisfactory"  link  flows 
{F  | aeL},  and  at  the  same  time  maintains  a  certain  degree  of  "fairness" 
for  all  sessions. 

It  is  customary  to  consider  as  one  of  the  characteristics  of  a  fair 
allocation  of  resources  in  a  network  the  feature  that  it  is  indifferent  to 
the  geographical  separation  of  the  session's  origin  and  destinations.  Al¬ 
though  there  might  be  different  priorities  assigned  to  sessions,  these 
priorities  are  not  assigned  on  the  basis  of  geographical  distance.  More¬ 
over,  two  sessions  of  the  same  priority  should  obtain  the  same  rate,  if 
the  rate  of  one  can  be  traded  for  the  rate  of  the  other  without  overload¬ 
ing  the  network  or  reducing  the  rate  of  any  other  session.  This  is  in 
the  spirit  of  making  the  network  "transparent"  to  the  user. 

To  capture  the  notion  of  fairness  and  priority  as  presented  above 
we  define  the  notion  of  fair  allocation: 
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For  a  vector  x  =  (x^.x2, . . .  .x11)  in  the  Euclidean  space  Rn,  we  consider 
the  vector  x  =  (x*,x2, . . .  .x11)  the  coordinates  of  which  are  the  same  as 
those  of  x  but  are  rearranged  in  order  of  increasing  value,  i.e.  we  have 
x1  <  x2...^x  nand  with  each  i  -  i,...,n  we  can  associate  a  distinct  i'  such 

_ i  i»  _ 

that  x  =«  x  .  We  call  x  the  increasing  permutation  of  x.  Given  a  subset 
X  of  Rn  we  will  say  that  a  vector  xeX  is  a  fair  allocation  over  X  if  for 
every  vector  yeX  the  increasing  permutation  x  of  x  is  lexicographically 
greater  or  equal  to  the  increasing  permutation  y  of  y,  i.e.  if  yJ  >  x? 
for  some  j,  then  there  exists  an  i  <  j  such  that  y1  <  x1. 

If  we  view  X  as  a  "feasible"  set,  a  fair  allocation  vector  x  over  X 
solves  a  hierarchy  of  problems.  The  first  problem  is  to  maximize  the 
minimal  coordinate  of  vectors  in  X.  The  second  problem  is  to  maximize  the 
second  minimal  coordinate  over  all  vectors  which  solve  the  first  problem, 
etc. 

Hayden  [2]  proposed  an  algorithm  which  results  in  a  rate  vector 
r  =  (...,  rg,  ...)  which  is  a  fair  allocation  over  the  set  defined  by 


Fa  <  p  Ca,  V  aeL,  (2) 

where  p  is  some  constant  between  0  and  1 .  Jaffe  £3]  proposed  an  algorithm 
which  obtains  a  rate  vector  r  such  that  the  vector  (...,  6srs*  •••)  is  a 
fair  allocation  over  the  set  defined  by 


6  r  <  C  -  F  ,  V  seS,  aeL  , 

s  s  —  a  a  '  *  s 


F 

a 


-  a' 


rg  >  0,  V  se$. 


(3) 


U 


t 


V  aeL 


(4) 
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where  gg  is  some  positive  constant  that  characterizes  the  priority  of 
session  s. 

The  rationale  behind  the  fair  allocation  problem  based  on  (2)  is  quite 
simple:  we  maximize  the  minimum  session  rate  while  not  allowing  the  flow 
of  any  link  to  be  more  than  some  given  fraction  of  its  capacity.  The 
rationale  behind  (3),  (4)  is  somewhat  more  sophisticated.  Primarily  it 
enables  us  to  establish  preferences  among  sessions,  and  to  accomodate 
fluctuations  of  a  session  rate  which  depend  linearly  on  the  rate  as  we  will 
demonstrate  shortly. 

While  Jaffe's  algorithm  is  not  iterative  and  as  a  result  is  somewhat 
unsuitable  for  distributed  operation,  Hayden’s  algorithm  may  result  in 
transient  flows  that  are  larger  than  some  link  capacities  (for  an  example 
see  [4],  p.  39). 

Our  purpose  in  this  paper  is  to  propose  and  analyze  an  iterative 
algorithm  that  solves  a  problem  that  is  more  general  than  Jaffe's  [3], 
maintains  at  all  times  feasibility  of  link  flows  with  respect  to  capacities, 
and  is  suitable  for  distributed  operation.  To  this  end  we  generalize 
the  set  defined  by  (3),  (4)  as  follows: 

For  each  link  aeL  and  session  seS  let  g  :  R+  R+  and  8  :  R+  -+•  R+ 

a  s 

be  functions  mapping  the  nonnegative  portion  of  the  real  line  R+  into  it¬ 
self.  We  are  interested  in  finding  a  rate  vector  r  such  that  the  vector 
(• • • *8s(rs) , • • •)  is  a  fair  allocation  over  the  set  defined  by 

W  -  v  seS’  aeLs  (5> 

F  <  C  ,  V  aeL,  r  >  0,  V  seS.  (6) 

a  “*  a  5  *“ 
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A  vector  r  with  this  property  will  be  called  a  fair  allocation  rate. 

We  make  the  following  assumptions  regarding  the  functions  g  (•)  and 

ci 

esc-): 

Assumption  A:  For  all  a eL ,  g  (•)  is  monotonically  nondecreasing,  and, 
for  all  seS,  8S(*)  is  continuous,  monotonically  increasing,  and  maps  R+ 
onto  R+.  (This  implies  also  that  the  inverse  g”1^)  exists,  is  con¬ 
tinuous,  monotonically  increasing  and  maps  R+  onto  R+) . 

Assumption  B:  The  function  Hsa(*)  defined  by 


Hsaco  -  s;1  [*,(«]. 


V  seS,  a eL,  feR 


is  convex  and  differentiable  on  R  and  satisfies 


Hsa(°)  =  0. 


Assumption  B  is  not  very  restrictive.  It  is  satisfied  in  particular 
if  hoth  A  (•)  and  g  (•)  are  convex,  differentiable  and  monotonically  in- 

S  a 

creasing  on  R+,  and  g  (0)  =  0.  Also  the  convexity  assumption  in  Assumption  B  can  be 
replaced  by  a  concavity  assumption  without  affecting  the  convergence 
result  of  the  next  section,  but  this  will  not  be  pursued  further. 

The  introduction  of  the  nonlinear  function  8S(0  allows  us  to  assign 
different  priorities  to  different  sessions  in  a  more  flexible  manner  than 
in  (3),  and  allows  additional  freedom  in  mathematically  expressing  algo¬ 
rithmic  design  objectives.  As  an  example  let  us  provide  justification 
for  the  use  of  a  particular  form  for  g  in  the  case  where  each  session 
is  a  voice  conversation. 
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Suppose  that  the  length  of  time  over  which  each  session  rate  is 
averaged  is  short  relative  to  the  "time  constant"  of  the  counting  process 
of  the  number  of  off-hook  speakers  which  are  currently  in  talkspurt  mode. 
Since  about  30%  of  a  talkspurt  is  silence  and  some  segments  of  the  talk- 
spurt  need  more  encoding  than  others,  we  view  the  bit  rate  generated  by 
the  Vocoder  for  session  s&S  as  a  stochastic  process  with  mean  rg — the  rate 
assigned  to  user  s&S.  We  thus  implicitly  assume  that  the  Vocoder  has  the 
means  of  dynamically  reconfiguring  to  the  demands  of  the  voice  to  achieve 
the  desired  average  rate.  Suppose  that  we  want  to  reserve  excess  capacity 
on  each  link  so  as  to  be  able  to  accomodate  a  variation  at  least  as  large 
as  the  standard  deviation  of  the  flow  on  the  link.  Assume  that  the 
standard  deviation  of  the  rate  of  each  session  s&S  is  y r  where 
0  <  y< 1 •  For  a  fixed  link  a eL  let  s'eS  be  such  that 

s'  =  arg  max  r  . 

s&Sa  S 
a 

Then,  by  the  independence  of  the  rates  of  different  sessions,  we  have, 
assuming  F  <  C  ,  that  the  standard  deviation  a(F  )  of  the  flow  F 
satisfies 


I 
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»»■  trs'>  ■  's'  •  *.(<VFa>  ■  -TT  ‘W 

Y  C. 


Then  from  (5),  (7)  and  (8)  we  obtain 


a(Fa)  <  Y^Cars,  <  Y^WW  “  Ca  ’  Fa' 


We  are  thus  guaranteed  to  be  able  to  accomodate  the  standard  deviation 
of  the  flow  resulting  from  the  fair  allocation. 

In  the  second  interpretation,  the  length  of  time  over  which  the 
rate  is  averaged  is  relatively  long  with  respect  to  the  "time  constant" 
of  the  counting  process  of  the  number  of  off-hook  speakers  in  talkspurt 
mode.  In  this  case  we  deal  concurrently  with  all  the  off-hook  sessions 
and  want  to  be  able  to  accommodate  the  standard  deviation  around  the 
mean  of  the  process  (i.e.,  the  instantaneous  effect  of  the  number  of 
speakers  at  the  talkspurt  mode  is  washed  out  by  the  long  time  average) . 
Let  q  be  the  fraction  of  time  a  speaker  is  in  the  talkspurt  mode  and 
assume  his  rate  while  in  the  talkspurt  mode  is  constant.  Then  using 
notations  as  before 


0(FJ  =  [  l  (r  /q)2*qCl-q)]1/2 


<  C— q  C  r*)1/2 

—  q  a  s' 


1  ^V1/2  ua(ca-Fa)i1/2- 
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Again,  by  choosing  g&  as  in  (8)  with  y  =  we  obtain 

1  VV 

The  point  we  want  to  make  by  the  above  arguments  is  that  there  is 
often  a  need  to  allow  ga  to  be  a  nonlinear  function,  which  may  depend  also  on 
Ca’  rather  than  only  on  the  excess  capacity  as  (3)  implies.  The  exact 
role  of  ga  is  up  to  the  network  designer  to  decide,  and  our  formulation 
allows  him  a  great  deal  of  flexibility  in  this  regard. 

It  is  possible  to  show  that  Assumptions  A  and  B  guarantee  existence 
.iffliff»eness  of  a  fair  allocation  rate.  The  proof  given  below  is  con¬ 
structive  and  is  based  on  a  finitely  terminating  algorithm.  However  this 
algorithm,  in  contrast  with  the  one  of  the  next  section,  is  not  suitable 
for  distributed,  on-line  operation  since  it  must  be  restarted  each  time 
an  old  session  is  terminated  or  a  new  one  is  initiated. 

Consider  first  the  problem  of  finding  a  vector  r  =  (...,  r  ,...)  that 
maximizes 


min 

seS 


MV 


over  the  feasible  set 

R0  =  tr[  (5)  and  (6)  are  satisfied} 


This  is  the  first  problem  in  the  hierarchy  of  problems  solved  by  a  fair 

allocation  rate,  and  can  be  solved  simply  by  observing  that  its  optimal 

value  [i.e.  max  min  B_(rJ]  is  equal  to 
teRq  seS 

W1  =  max{w|w  <  ga[ca  -  J  B^OOJ,  ael}  . 


(9) 
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This  follows  easily  from  the  fact  that  both  g  and  g  *  are  monotonically 

3.  S 

nondecreasing.  Denote 

L*(  1)  -  {aei|w*  *  g  [C  -  J  glV?)]} 

1  a  a  s  eSa  S  1 

S*(l)  =  {se5  |  Lsm*CD  *  0). 

For  any  fair  allocation  rate  (...,  r  ,  ...)  the  rate  of  the  sessions  in 
S’* Cl)  is  equal  to  g^Cw*),  i.e. 

rs  =  VseS*(l), 


while  L* Cl)  may  be  viewed  as  the  set  of  bottleneck  links  the  presence 

of  which  does  not  allow  us  to  increase  min  g  (r  )  beyond  the  level  w* 

seS  S  5  1 

Therefore  for  the  purposes  of  determining  further  a  fair  allocation 
rate  vector,  the  rates  of  the  sessions  seS*(l)  are  fixed  at  g^(w*)  and 
we  can  consider  a  reduced  network  whereby  the  links  aeL*(l)  and  sessions 
seS*(l)  are  eliminated  while  the  capacity  C  of  each  link  a£L*(.l)  is 
replaced  by 


C 

a 


sLnnns. 


If  S* Cl)  *  S  we  are  done;  otherwise  we  can  consider  the  problem  of 
maximizing  the  minimal  value  of  Ss(rs)  in  the  reduced  network  similarly 
as  earlier.  This  will  determine  a  new  optimal  value  w*  with  w*  >  w*,  a 
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a  new  set  of  bottleneck  links  L*C 2),  and  a  set  of  sessions  S*(2)  such  that 
Ss(rs)  =  «*  ,  V  seS*(2) 

in  any  fair  allocation  vector.  If  S*(1)US*(2)  =  S  we  are  done;  otherwise 
we  can  proceed  by  constructing  a  reduced  network  and  continue  in  the  same 
manner  as  earlier  until  we  exhaust  all  sessions.  This  argument  constructs 
a  fair  allocation  rate  r*  and  shows  that  it  is  uniquely  defined  in  terms  of 
the  scalars  wj,  w*,...,  and  the  corresponding  sets  S*(l)>  S*(2),...  Note 
also  that  the  session  rates  r*»  s&S  and  associated  link  flows  F*, aeL 
satisfies 

r*  =  min  H  (C  -F*) ,  V  se S-  (10) 

s  a CL  sa  a  a 

s 

The  algorithm  of  the  next  section  is  based  on  this  property. 

We  can  also  show  the  reverse  property  namely  that  if  a  rate  vector 
r*  satisfies  (10)  then  it  is  a  fair  allocation  rate.  To  see  this  let 
r  =  (...,r  ,...)  satisfy  (10)  or  equivalently 

6s(rs)  =  min  8a(Ca-Fa),  V  seS.  (11) 

aeL 

s 

Observe  that  from  the  definition  (9)  of  wj  and  (11)  we  obtain 


min  min  g_(C  -F). 
seS  aeL 

s 


where 


Let  a eL  be  any  link  such  that 


g~ (C~-F~)  =  w, . 

sa  a  a  1 


Then  from  (11)  we  have 


?s  =  *s  ^l3’  V  seSa* 


The  inequality  w.  <  w*  implies  that  F~  =  ?  r  <  £  r*  =  F£ 

a  sES~  s  "  seS~  s  a 

where  r*  is  the  fair  allocation  rate.  But  this  implies  that 


*i  =  *s<crFs>  i  H«hr*V  >  WI- 


Therefore  we  must  have  w^  =  w*  and  it  follows  that  the  vector  ?  solves 
the  first  problem  in  the  hierarchy  of  the  fair  allocation  problem. 
Proceeding  similarly  as  earlier  we  can  show  that  r  solves  all  the  problems 
in  the  hierarchy  of  the  fair  allocation  problem  and  is  therefore  a  fair 
allocation  rate. 

We  summarize  the  conclusions  from  the  preceding  arguments  in  the 
following  proposition. 

Proposition  1;  Let  Assumptions  A  and  B  hold. 

a)  There  exists  a  unique  fair  allocation  rate. 

b)  r*  »  (...,r*,...)  is  a  fair  allocation  rate  if  and  only  if  it 


satisfies 


Bs(r*)  =  min  ga(Ca-F*),  V  seS 
ael 

s 


or  equivalents 
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T1  =  “in  V  S£S- 
a cLs  sa  a  a 


(13) 


In  some  situations  it  may  be  reasonable  to  consider,  in  addition  to 
(5)  and  (6) ,  the  constraint 


rs  £  Rs.  V  seS 


(14) 


where  Rg  is  given  upper  bound  to  the  rate  of  session  rate  s.  We  may  view 

Rs  either  as  a  limit  on  rate  imposed  by  technological  restrictions  or 

as  a  maximum  desired  rate  by  session  s.  The  problem  of  finding  a  fair 

allocation  rate  over  the  set  defined  by  (5),  (6),  and  (14)  can  be  reduced 

to  the  problem  considered  earlier  by  introducing  for  each  seS,  an 

artificial  link  ag  traversed  only  by  session  s  by  setting  the  capacity 

Ca  of  that  link  equal  to 
s 

%  ■  RS*W- 

and  by  selecting  the  function  g  to  be  the  identity.  Then  the  constraint 


Mrs>  i  «a  (ca  -rJ 
s  s 


S'- s'  -  °a_"X  ‘s'  -  ca  -  rs 


becomes  rg  +  8g(rg)  <  *  Rg  ♦  8g(Rg)  ar*d  Is  equivalent  to  (14). 

s 
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3.  THE  ALGORITHM 

k  k 

Let  r  =(...,  rg,  . ..)  be  the  rate  vector  obtained  after  k  iterations 

k 

and  let  {F  }  be  the  corresponding  set  of  total  link  flows.  Assume  that 

SL 

0  <  F J  <  C.  V  ael.  (15) 

a  a 


k+l  k+l 

The  new  rate  vector  r  *(•••»  r  *••■)  obtained  at  the  (k+l)st  iteration 
is  given  by 


k+l 


mm 

aeU 


k  k  k 

{r  +  v  [H  (C  -Fk) 
1  s  ’a  1  sa  a  a 


ts)},  V  seS 


(16) 


where  ya  is  given  by 


k 


1  +  y  H*  (C  -Fk  ) 

s is  sa  a  a 


(17) 


and  H'  (•)  denotes  the  first  derivative  of  H  (•)•  In  a  practical 
sa  sa 

k 

implementation  of  the  algorithm  the  link  flows  F  can  either  be  measured 

a 

(by  taking  a  time  average),  or  they  can  be  mathematically  computed  as 

k 

the  sum  of  the  session  rates  r^,  s&Sa-  The  session  rates  computed  via 

(16),  (17)  will  have  to  be  translated  into  physical  rates  by  software 

residing  at  the  session  origin  nodes. 

The  following  lemma  shows  among  other  things  that  property  (15)  is 

maintained  by  the  algorithm  at  each  iteration  and  therefore  if  the  initial 

link  flows  are  within  the  link  capacities  C  the  same  is  true  for  all 

& _ 

link  flows  generated  by  the  algorithm. 

Lemma  1 :  Let  Assumptions  A  and  B  hold  and  assume  that  the  initial  rate 
vector  r°  is  such  that 
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i 

i 

i 

i 


0  <  r°  ,  *  Ses,  0  <  F!  <  C„.  V  a el.  (18) 

5  ft  ft 

k  k 

Then  if  {r  }  is  a  rate  sequence  generated  by  the  algorithm  (16)  with  y 

ft 

given  by  (17)  we  have  for  all  k 

0  <  r*  ,  V  seS,  0  <  Fk  <  C&,  V  aei.  (19) 

Furthermore 

Fk  <  7  H  (C  -Fk),  V  a eL,  k  >  1.  (20) 

a  —  £  sa  a  a  — 

seS. 


Proof :  See  Appendix  A. 


The  idea  behind  the  choice  of  expression  (17)  as  well  as  the  intuition 
behind  Lemma  1  can  be  best  explained  by  the  use  of  Figure  1.  Let  the  function 
G  (•):  R+  ■*  R+  be  given  by 

ft 


W  *  l.  Hsa<Ca-V 


seS. 


sa'  a  a' 


(21) 


k  1 

The  figure  depicts  the  relations  between  F a  and  Fa  ,  as  if  the  network 

k+1 

consisted  of  the  single  link  a.  F  is  determined  by  intersection  the 

ft 

k  k  •  • 

tangent  to  the  graph  of  G&(Fa)  at  the  point  (Fa>  with  the  line 


k 

y  =  Fa.  The  reader  can  easily  convince  himself  that  lim  sup  Fa  must  lie 

k-K» 

in  the  area  where 


Fa  i  W 

which  gives  rise  to  the  lemma.  Figure  2  shows  why  just  monotonicity  of 
G  (•)  is  not  sufficient  for  the  lemma  to  hold. 

ft 

We  can  now  state  the  main  result  of  this  paper. 


If 

Proposition  2:  Under  the  assumptions  of  Lemma  1  the  sequence  {r  }  con¬ 
verges  to  the  fair  allocation  rate. 


Proof :  See  Appendix  A. 

A  Variation  of  the  Algorithm. 

In  iteration  (16)  we  have  assumed  that  updating  of  all  the  rates 
rs  takes  place  simultaneously.  It  is  possible  to  consider  other  related 
algorithms  whereby  a  single  session  rate  rg  is  updated  using  (16) ,  then 
the  flows  F  are  updated  to  reflect  the  change  in  r  ,  then  another  session 
rate  is  updated  using  (16)  and  so  on  until  all  the  session  rates  are 
taken  up  cyclically  in  a  fixed  order.  This  one-session-at-a-time  mode 
of  operation  is  reminiscent  of  the  Gauss-Seidel  method  for  solving  systems 
of  equations  and  is  perhaps  better  suited  for  distributed  implementation. 

It  is  possible  to  show  that  all  the  results  of  this  section  hold  for  this 
modified  algorithm  as  well. 

A  question  of  considerable  interest  is  whether  a  totally  asynchronous 
distributed  version  of  algorithm  (16)  will  work  satisfactorily  (compare  with 
algorithms  investigated  in  [5],  [6]).  In  such  an  algorithm  each  session  orivin 
sends  at  arbitrary  times  along  the  session  path  a  control  packet  contain¬ 
ing  the  current  rate  of  the  session.  As  the  packet  travels  to  its 
destination  the  information  needed  to  compute  the  right  side  of  (16)  is 
collected.  (We  assume  here  a  single  destination  per  session  and  that 

each  link  a  on  the  session  path  maintains  the  current  value  of  F  as  the 

a 

sum  of  all  currently  assigned  session  Tates  r  ,  seS  ,  and 

s  a 
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the  fora  of  the  function  H  (•)  for  each  seS  ) .  The  destination  returns 

SE  cl 

the  new  rate  to  the  session  origin  and  the  links  along  the  session  path. 
This  type  of  algorithm  is  very  attractive  from  the  practical  point  of 
view  since  it  does  not  require  a  session  synchronization  protocol.  Its 
convergence  properties  however  are  as  yet  unclear  and  are  currently 
under  investigation.  It  is  interesting  to  note  that  some  of  the  algo¬ 
rithms  investigated  by  simulation  in  [1]  are  of  similar  nature. 
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APPENDIX  A 

Proof  of  Lemna  1; 

It  suffices  to  show  (19)  -(20)  for  k  =  1.  Consider  the  function 

G  :  R+  ■+  R+  defined  by 
a 


W 


*  l  H  (C  -FJ  if  0  <  F  <  C 
sav  a  a  —  a  — 


seS_ 


a  —  a 


if  Fa  >  Ca* 


From  (16)  and  (17)  we  have 


1  \  o 

r  =  min  It  +  - 

•cts  I  l-Ck(F°) 


lHaa'Ca-Fa>  '  rs>  |  '  T 


or 


min 

aeL, 


(Hsa<Ca-tf  -  \ 

\  i-g;cp °)  )  ' 


v  seS . 


(A.l) 


Since  G  (•)  is  monotonically  nonincreasing  we  have  G’(F°)  <  0,  and  since 
also  H  (C  -F°)  >  0  we  obtain  from  the  hypothesis  r°  >  0  and  (A.l) 

S3  3  <*  “ ’  S  “ 


l 

r 

s 


>  0, 


and  therefore  also 


V  seS, 


(A.  2) 


F 


1 

a 


>  0, 


V  aeL 


(A.  3) 
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From  (A. 1)  we  have 


H  (C  -F  )  -  r°G’(F°) 

co  o  n  e  e  *•  e' 


r  < 
s  — 


i  -  g;(fJ 


V  seS  ,  a eL 


and  by  adding  over  all  seSa  we  obtain 

.  G  (F°)  -  F°G'(F°) 

_1  .  a  a'  a  a  aJ 

F  <  - - -  , 

3  i  rti 


V  a  eL. 


1  ' 


Since  1-G'(F  )  >  0  we  obtain  from  the  inequality  above 

&  3 


F1  <  G  (F°)  +  (F1-F°)G' (F°) 
a  —  a.K  a.J  ^  a  a  av  a 


Since  Ga(")  is  convex  the  right  side  of  this  inequality  is  less  or  equal 
to  G  (F*)  and  we  obtain 

A  A 

<a  i  GaCFa>'  CA-4J 


Since  G  (•)  is  monotonically  nonincreasing  and  G  (F)  *  0  for  F  >  C  we 
a  a,  —  a 

obtain  from  (A. 4) 


F 


1 

a 


<  C  . 


a 


(A. 5) 


From  (A.2)-(A.5)  we  obtain  (19)  and  (20).  Q.E.D. 

Proof  of  Proposition  2; 

Denote 

r*  =  lim  inf  r^  ,  V  seS. 
s  k—  5 
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k 

Fix  s&S  and  consider  a  subsequence  {r  }.  v  converging  to  r*.  We  have 

S  S 

from  (16) 


r* 

s 


lim  min  {rk 

k-*»  a£/  s  a 

keKs 


[H  (C  -Fk_1) 
1  sav  a  a  1 


]}  • 


Since  L  consists  of  a  finite  number  of  links  we  may  assume  (by  passing 
to  a  subsequence  of  K  if  necessary)  that  there  exists  a  link  a  such 

S  5 


that 


*  r  k-1  ^  k-1  ru  .  ck-l  k-1, . 

r*  =  lim  {r  +  y  [H  (C  -F  )  -  r  ]} 

s  k-K  s  as  sas  as  as  s 
s 

k-x» 


(A.  6) 


k-1 

Since  {F  }  is  bounded  above  and  below  we  may  assume  (by  passing  to  a 
as  keK 

s 

subsequence  of  K  if  necessary)  that  for  some  F 

s  as 

k-1  " 

lim  FK  1  =  F  . 

■%  i/  a  a 

k-Ks  s  s 

k-HXJ 


Denote  also 


^  <Fa  ^ 
s  s 


=  lim  Y, 

ke<s 

k-*» 


k-1 


We  have  from  (A. 6) 


r*  >  ( 1 — Y  )  lim  inf  rk +  v  H  (C  -F  ) 

s  -  v  'a  s  'a  sa. v  a.  a  ' 


s  k-*-K 

k-«o 


s  s  s  s 


>  Cl-Y_  )r*  ♦  Ya  H  (C  -  lim  sup  Fy  ) 
as  s  s  s  s  s  k-K»  as 
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and  finally 


r* 

s 


>  H  (C  -  lim  sup 
s  as  k-*» 


(A.  7) 


Since  the  choice  of  s  was  arbitrary  we  conclude  that  for  every  s&S  there 

exists  a  eL  such  that  (A .7)  holds, 
s  s 

Let  a jet  be  such  that 


a,  =  arg  min  g  (C  -lim  sup  F*) 
i  «  a  a  •  a. 

ae  L  k-*» 


Using  the  monotonicity  of  g^1  we  obtain 


H  (C  -lim  sup  F  )  _>  H  (C  -lim  sup  F  ) 


sa  a  ,  *  a 

s  s  k-x»  s 


and  therefore  from  (A. 7) 


sa,  a,  ,  r  a, 

1  1  k-«o  l 


>  H  (C  -lim  sup  F*  ) ,  V  seS 


SV  al 


k-KO 


Summing  (A. 8)  over  all  seS  we  obtain 

al 


lim  inf  F  £  lim  inf 

k-x»  al  seS 

a. 


seS 


-  ..L  »  (C  -lim  sup  F*  ) 


ShTo  sa.  a 


1  “1  k-~» 


(A.  8) 


On  the  other  hand  from  (20]  we  have 
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lim  sup  Fk  <  I  H  (C  -  lira  sup  ). 
k  al  “  seS  sal  al  k  al 

al 


It  follows  that  the  last  two  inequalities  as  well  as  (A. 8)  hold  as  equations, 

the  entire  sequence  {F^  }  converges  to  £  r*  while  each  sequence 

al  seS  5 

k  &1 

{r  },  seS  ,  converges  to  r*. 
s  s 

Consider  now  a  new  network  derived  from  the  previous  one  by  deleting 
link  a^,  and  all  the  sessions  traversing  it.  We  consider  the  algorithm 
executed  in  the  same  manner  as  before  with  the  same  initial  rates  for 
the  remaining  sessions  but  with  the  capacity  of  each  link  ael  replaced  by 


-  I 

seS, 


This  will  result  in  the  same  rate  sequence  for  the  sessions  s£  S  as  in 

al 

the  original  algorithm.  A  trivial  modification  of  the  argument  used  to 
show  (20)  in  Lenma  1  shows  that  we  will  have 


lim  sup  F^  <  £  H  (C  '  lim  su^  v  ael  • 

k^»  a  ”  seS  sa  3  k+~  a 

a 


This  relation  can  be  used  to  repeat  the  argument  given  earlier  in  order 

k 

to  show  the  convergence  of  the  sequences  {rs>  to  r*  for  all  sessions 

s£  S_  traversing  the  link  a-  where 
al  i 

Ir 

a~  *  arg  min  g fC  -lim  sup  F  )  . 

1  aeL  a  a  k-H»  a 

a^ax 


By  repeating  this  procedure  we  will  eventually  exhaust  all  links 

y 

thereby  showing  that  each  rate  sequence  {r  },  seS  converges  to  r*,  each 

s  s 
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k 

each  flow  sequence  {F  },  ae L  converges  to  F*  =  7  r*  and  each  stepsize 

a  a  se$  s 

a 


sequence  {y  },  aeL  converges  to 


Y* 

'a 


1  *  l  H-<c^> 

a 


By  taking  limits 


r*  = 
s 


in  (16)  we  have  for  all  seS 


lira  min  {r*  +  y*[H  (C  -F*)  -  r*] } 
k-*»  aeL  S  a  sa  a  a  s 
s 


min  lim  {r^  +  y^[H  (C  -F^)  -  r  ]) 
.  .  s  TaL  sav  a  a  sJ 

a  cL  k-xo 
s 


(A.  9) 


where  the  interchange  of  min  and  lim  above  is  valid  since  all  the  sequences 
inside  the  braces  converge  and  the  number  of  elements  of  Lg  is  finite. 

From  (A. 9)  we  obtain  for  all  seS 


min 

aeL 


s 


^Hsa^FP 


r*3 


o 


Since  y*  >  0  for  all  aeL  we  obtain 
a  s 


ra  =  mi"  Hsa(Ca-Fa)  ’  V  seS< 
aeLs 


The  result  now  follows  from  Proposition  1  [cf.  (13)] 


Q.E.D. 
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